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y—{ Abstract 

^ In order to establish links between p-wave pion production in nucleon-nucleon collisions and low energy three-nucleon scattering, an 
k>( extensive programme of experiments on pion production is currently underway at COSY-ANKE. The final proton pair is measured 
at very low excitation energy, leading to an 5*- wave diproton, denoted here as {ppjs . By using a deuterium target we have obtained 

d data on the differential cross section and analysing power of the quasi-free pn — >■ {pp}s7r~ reaction at 353 MeV. The spectator 
proton psp was either measured directly in silicon tracking telescopes or reconstructed using the momentum of a detected vr" . Both 
observables can be described in terms of s-, p-, and d-wave pion production amplitudes. Taken together with the analogous data 
on the pp — >■ {pp}a tt" reaction, full partial wave decompositions of both processes were carried out. 
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There is an extensive programme of near-threshold mea- 
surements of NN {pp}s IT at the COSY-ANKE facility of 
the Forschungszentrum Jiilich [1,2]. Here the {pp\s denotes 
a proton-proton system with very low excitation energy, 
Epp, which is overwhelmingly in the ^5*0 state with antipar- 
allel proton spins. The primary aim of these experiments is 
to carry out a full amplitude analysis which would lead to a 
determination of the pion p-wave production strength from 
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the ^Si initial state that could provide links with other in- 
termediate energy phenomena. As part of this programme, 
we have already presented data on the cross section and 
analysing power of the pp — > {pp}s reaction, which can 
be described in terms of s- and d-wave pion production 
ampUtudes [3]. We here present analogous results for the 
quasi-free pn —> {pp}s7r~ reaction at an effective beam en- 
ergy around 353 MeV. 

In the absence of an intense monochromatic neutron 
beam, the study of the pn — >■ {pp}s7r~ reaction is most eas- 
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ily carried out through quasi-free 7r~ production, with a 

proton beam incident on a deuterium target. In order to 
make a full reconstruction of a — )■ Psp{pp}s tt" event, it 
is necessary to measure accurately the momentum of three 
of the particles in the final state. Determinations have been 
made at TRIUMF of both the cross section [4] and pro- 
ton analysing power [5] by detecting the pion together with 
the two fast protons. The spectator proton, psp, was then 
identified through the missing mass in the reaction and its 
momentum reconstructed kinematically. Data were taken 
for three beam energies, though only in the central region 
of pion angles. The lowest of these energies was 353 MeV 
and, to allow a close comparison, we have concentrated our 
attention on this energy. 

In contrast to the TRIUMF data, we have extended the 
coverage to the whole angular domain, which is very impor- 
tant in the subsequent partial wave analysis. Furthermore, 
by including the possibility of detecting directly the spec- 
tator proton or the 7r~ , we have obtained an independent 
check on some of the systematics involved. 
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Fig. 1. Top view of the ANKE spectrometer setup, showing the 
positions of the Positive (PD), Negative (ND), and Forward (FD) 
detectors, as well as the Silicon Tracking Telescope (STT). The 
dipoles Dl and D3 deflect the circulating proton beam in and out 
of ANKE, whereas D2 serves as an analysing magnet. 

The experiment was carried out at the ANKE spectrom- 
eter facility [6] , which is installed inside the COSY cooler 
synchrotron storage ring of the Forschungszentrum Jiilich. 
The circulating proton beam was polarised perpendicularly 
to the horizontal plane of the machine and the polarisation 
reversed in direction every six minutes. 

Fast protons arising from the interaction of the beam 
with the deuterium cluster-jet target [7] traversed the spec- 
trometer dipole magnet D2 shown in Fig. 1 and entered the 
forward (FD) and/or positive side (PD) detectors. Negative 
pions produced in the interaction at small centre-of-mass 
(CM) angles, 6*^ < 40°, could be detected in the negative 
side detector (ND), while slow spectator protons (psp) were 
observed in one of the silicon tracking telescopes (STT) 



that were located in the target vacuum chamber close to 

the target jet [8,9]. 

Each of the FD, PD, and ND detectors shown in Fig. 1 
includes both scintillation counters and multiwire propor- 
tional (MWPC) or drift chambers (MWDC). The recon- 
struction of the particle trajectories from hits in the wire 
chambers allowed the momenta of the ejectiles to be evalu- 
ated. The counters were used to measure the arrival times 
and energy losses required for particle identification. 

The two STT were placed symmetrically to the left and 
right of the target jet. Each telescope consists of three sen- 
sitive silicon layers, though only the first two were used 
in this experiment. Spectator protons with the energies 
2.5 < -Esp < 6 MeV pass through the first layer and 
stop in the second. The accuracy of the energy measure- 
ment, a{Esp) / Esp — 10%, allow protons to be identified by 
their energy loss [10] . Although higher energy recoil protons 
could be measured by using the third STT layer, to stay 
within the range of applicability of the spectator model, 
only events with Egp < 6 MeV were retained. It should, 
however, be noted that there is no lower limit on Egp when 
the 7r~ is measured. 

In order to provide suflacient resolution in Epp, both fast 
protons from the diproton have to be detected and, to per- 
mit a complete reconstrtiction of the kinematics of the re- 
action, one has in addition to detect either the 7r~ or the 
slow spectator proton. The triggers selected fast proton 
pairs, with either both protons hitting the FD or PD, or 
one proton recorded in each detector. A trigger for single 
tracks in the FD was also added to record deuterons from 
the pn — >■ dn^ process, which were used for normalisation 
and polarimetry purposes. A sample of data was also taken 
with a special trigger requiring a coincidence between the 
first two layers of the STT. This was used to investigate 
the STT performance and for background studies. 

The CM energy of the quasi-free pn system depends on 
the energy and the emission angle of the spectator and 
for this experiment the effective "free" beam energy was 
in the range Tfreo ~ (310-390) MeV. However, for the re- 
constructed events. Tficc was measured with an accuracy 
of a{T{ree) = 2-4 MeV and only data in the range Tfree = 
353 ± 20 MeV were used in the analysis. 

The identification of the pd ^ Psp{pp}sT^^ reaction 
starts from the selection of pairs of fast protons through 
the difference between their times of ffights in ANKE [11]. 
Since the beam energy is qtiite low, other pairs of par- 
ticles gave a negligible contribution to the background, 
which originated mostly from accidental coincidences, and 
amounted to < 7%. To select the; diproton state, we 
imposed a Epp < 3 MeV cut on the data, which could 
be done reliably because of the excellent resolution of 
a{Epp) < 0.6 MeV. It is important to note that the iden- 
tical cut was placed on the pp — ^ {pp}s'^° data [3], which 
is important when considering the relative normalisation 
uncertainties for 7r° and tt" production. 

The time-of-flight criterion was also used to select the tt" 
in the ND, where the very low background consisted only of 



accidentals and strongly scattered positively charged par- 
ticles. The background level for the spectator protons iden- 
tified in the STT was at the 5-8% level, depending on the 
proton energy and angle. 

Having measured three of the final particles, the resid- 
ual one from the pd — > Psp{pp}s'!T^ reaction was identified 
by the missing-mass method. This is illustrated separately 
in Fig. 2 for the cases where the spectator proton or 7r~ 
was detected. The main backgrounds were from accidental 
coincidences, which were particularly significant for spec- 
tator protons because no timing information from the STT 
was used in the analysis. The shape of this background 
was derived from artificially constructed events. The mo- 
menta of the fast protons were taken from the experimen- 
tal events where neither a spectator proton nor tt^ was de- 
tected, to which was added a random spectator momentum 
taken from the sample acquired with the STT trigger. The 
missing-mass spectrum was then fitted with the sum of this 
background distribution and a Gaussian. The background 
level was estimated separately for each detector combina- 
tion, spin orientation of the beam, and angular bin. Any un- 
certainties here were combined with the statistical errors. 
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Fig. 2. Kinematic identification of the pd — >■ Psp{pp}s 7r~ rea<;tion. a) 
Sample of missing-mass squared data from events where the spectator 
proton is detected, showing the experimental spectrum (with error 
bars), the background (shaded area), and the sum of this plus a 
Gaussian for the peak (solid curve), b) Missing-mass distribution 
when the is detected. The curve shows the fit to the experimental 
data with the sum of a Gaussian centred on the mass of the proton 
and a linear ba<;kground. 

In order to estimate the resolutions and acceptance of 
the system, a full simulation of the ANKE setup was per- 
formed [11], based on the GEANT package [12]. The ac- 
ceptance was calculated as a function of 9,^ for both the ■jt~ 
and spectator proton detection, assuming that the specta- 
tor was emitted isotropically, with the energy being sam- 
pled from the Fermi distribution [13]. The fast proton pair 
was generated in the ^Sq state, with the distribution in ex- 
citation energy being weighted by the Migdal- Watson fac- 
tor [14] that included the Coulomb interaction [15]. 

The polarisation of the proton beam and the luminosity 
were both estimated from quasi-free pn dn^ data that 
were taken in parallel. The fast deuteron was detected in the 



FD, and selected by its energy loss in the counters, and the 

spectator proton measured in coincidence in the STT. The 
reaction was then identified by the missing-mass method. 

The measurements of both the analysing power and cross 
section depend sensitively on the relative luminosities for 
the two spin states and this was controlled by comparing 
the rates of ejectiles emitted at 6* = 0° or = ±90° , which 
cannot depend upon the beam polarisation. Using calibra- 
tion data taken from the SAID database [16], the integrated 
luminosity was found to be L = (2310 ± 110) nb"''^. No cor- 
rection was made for the shadowing in the deuteron since 
there is a similar effect in the measurement of the quasi-free 
pn — ^ {pp}s'^~ rate. The beam polarisation determined in 
this way was P = 0.66 ± 0.06. This is consistent with that 
found ill the tt" production experiment [3] that was under- 
taken with the same conditions just before the 7r~ run. The 
luminosities in the two polarisation states were on average 
very close, L^/L; = 1.017 ± 0.005. 

Figure 3 shows the differential cross-section da/dil, for 
the quasi-free pn — > {pp}sTT~ reaction integrated over the 
Epp = (0 — 3) McV range oipp excitation energies and aver- 
aged over the effective beam energy Tfree = 353 ± 20 MeV. 
The data were extracted within the impulse approximation 
model, where the weighting of the spectator momentum 
distribution was taken from the Bonn deuteron wave func- 
tion [17], though the result was insensitive to this choice. 
The cross-section was estimated separately for tt" and spec- 
tator proton detection and, since these two results were 
found to be consistent, only their weighted average is pre- 
sented. Also shown are the TRIUMF data on quasi-free tt" 
production [4]. They imposed a slightly more severe Epp 
cut on their data and these were converted to our 3 MeV 
cut by assuming a Migdal- Watson energy variation [14]. 

Whereas the TRIUMF results only cover the; central re- 
gion of pion angles [4], the current data extend over the 
whole angular domain. The two data sets are consistent in 
the backward hemisphere but the TRIUMF measurements 
show no indication of the rise at forward angles that is seen 
at ANKE. Some confirmation of the ANKE angular shape 
is offered by pion absorption data, 7r~^He pnn^p, where 
the unobserved slow neutron is assumed to be a specta- 
tor [18] . In this case the reaction can be interpreted as being 
'!T~{pp}s — > pn, though the internal structure of the dipro- 
ton is very different to that in the production data. Over 
the range of angles covered, our data are completely consis- 
tent with these absorption results. The forward/backward 
peaking is in complete contrast to the results found for 7r° 
production [20,3] and is an indication of the dominance of 
the 7 = p-wave amplitudes in this reaction. 

The unpolarised cross section for tt~ production, and this 
times the proton analysing power Ay, must be of the form 

(da \ k 
^J^ = ^sine^^6„+icos"^„, (2) 
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Fig. 3. Unpolarised difFerential cross section for the pn — > {pp}sTT~ 
reaction at Tfroc ~ 353 MeV. The ANKE data with statistical errors 
are shown by red circles. In addition there is a systematic error of 
6% arising from the luminosity and acceptance determination. The 
statistical errors of the TRIUMF pn — > {pp}s7r~ results [4] (green 
triangles) are generally smaller than the symbol size and their nor- 
malisation uncertainty is 10%. The arbitrarily scaled TRIUMF cross 
sections extracted from 7r~^He — > pnnsp data [18] (blue stars) are 
also included. The dashed curve is a direct cubic fit to these ANKE 
data whereas the solid one corresponds to the global fit described in 
the text. 

where is the pion cm. production angle with respect 
to the direction of the polarised proton beam. Here p is 
the incident cm. momentum and k that of the produced 
pion. We are neglecting here any small effects due to the 
mass differences and at 353 MeV; the momenta then have 
values p = 407 MeV/c and fc « 94 MeV/c. The best fits of 
Eq. (1) to the differential cross section are found with the 
parameters quoted in Table 1. 



Observable 


Direct fit 


Global fit 


ao(PP) 


4.05 ±0.08 


4.05 ±0.08 


02 (pp) 


-2.31 ±0.14 


-2.34 ±0.14 


62 (pp) 


1.82 ±0.10 


1.80 ±0.10 


ao(pn) 


2.69 ±0.18 


2.47 ±0.08 


ai (pn) 


-8.24 ±0.51 


-7.83 ±0.45 


12 (pn) 


9.11 ±0.70 


10.12 ±0.41 


a-iipn) 


2.89 ±0.90 


1.38 ±0.27 


bi(pn) 


1.77 ±0.14 


1.82 ±0.13 


62 (pn) 


-1.95 ±0.50 


-1.75 ±0.36 


ba (pn) 


-4.43 ±0.70 


-4.83 ±0.27 



Table 1 

Values of the parameters in fih/sr extracted by direct fits of Eqs. (1) 
and (2) to the pp — > {pp}s tt" and pn — >■ {pp}s7r~ experimental data 
and those obtained on the basis of Eq. (5) from the amplitudes given 
in Eq. (6). The error bars are purely statistical. In the 7r~ case there 
are systematic uncertainties of 6% in the cross section and 9% in the 
analysing powers. 

The results for the analysing power of the pn {pp}s7r^ 



reaction are displayed in Fig. 4, with Ay{d(j/dil) being 
shown in panel a and Ay in panel b. The agreement with 
the TRIUMF Ay data [5] is reasonable at large angles and 
both show the strong and rather asymmetric oscillation in 
the central region. However, there are clear discrepancies 
for 9-^ < 60°, as there are also for the cross section shown 
in Fig. 3. 
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Fig. 4. (a) The product of the measured analysing power and dif- 
ferential cross section for the pn — > {pp}s7T~ reaction at 353 MeV, 
where the error bars shown are statistical and do not include the 11% 
systematic uncertainty. The dashed curve represents the best fit of 
Eq. (2) with hi, b2, and 63 terms whereas the solid one corresponds 
to the global fit to all the data, (b) Measured values of Ay for the 
pn — >■ {pp}s7T~ reaction showing both the ANKE (circles) and TRI- 
UMF data [5] (triangles). The systematic uncertainty in the ANKE 
data is 9%. The lines represent the quotients of the fits in panel-a 
and those to the cross section in Fig 3. 

Fitting the weighted Ay distribution with the form of 
Eq. (2) requires at least the three terms given in Table 1. 
The associated curve is shown in Fig. 4a and this divided 
by the parameterisation of the cross section in panel b. 

The spin structure of the pn — {pp}sTT~ reaction is iden- 
tical to that of pp — >■ {pp}s 7r° , which was already presented 
in [3] . The cross section and analysing power can be written 
in terms of two scalar amplitudes A and B through 



A,. 



da 

dn 

da 

dn 



Ap 



{\A\^ + \B\'^ + 2Re[AB*]cose^), 



= — (2ImUB*]sin6'„) 
4p 



(3) 



Keeping terms up to pion d waves, the pp — )■ {pp}s tt'^ 
data at 353 MeV [3] can be parameterised in terms of the 
three amplitudes , , and , corresponding to the 
transitions, ^Pq ^Sqs, ^Sod, and ^i^2 -> ^Sod, 



A 



respectively. In proton-neutron collisions there are also the 
two p-wave transitions, ^Si — > ^Sqp and ^Di — > ^Sqp that 
arise in the isospin / = case, and for these we introduce 
amplitudes and Mp, respectively. 

In terms of these partial waves, A and B become 

A=^ [Mf - iMj- + Mi (cos^ e^-l)+M^ cos 6^] , 

B = ^ [M^ - + {MP - IM!) cos e^] . (4) 

It is clear from the very different behaviour of the angular 
distributions for 7r° and 7r~ production that the extra p- 
wave amplitudes in Eq. (4) must be very large and that it 
will not be justified to discard their interference with M^ . 

Keeping terms up to p-d interference but omitting the 
squares of the c?-wave amplitudes, we find that 

ao = \\M^\^ + i|Mp^ - iMp^p - iRe [Mf*iM[ + f M,^)] 
ai = Re [Mf *(Mp« + ^M^) + ^MtiM^ - ^M^) 
-iMtM^] 

a2 = l\M^\' + Re [Mf *(Mf + |M,^) + M^*M^] 
as = Re [[M^*{Mi + ±M[) + M^*M[] 
h=lm [{M^* - \M^*){Ms - \{Mi + |M<f ))] 
62 = Im [M^*M^ - M^*{Mr, U4)] 
bs=lm[MtM^ - Mt{M^ + ^M^)] . (5) 

By neglecting the small coupling between the and 
partial waves, and imposing the Watson theorem to deter- 
mine the phases, it was possible to extract values for the 
complex amplitudes Mf, , and Mj" from the analysis 
of the pp — >■ {pp}s 7r° reaction [3] . Such an approach would 
not be valid for the two p-wave terms because of the very 
strong coupling between the incident ^Si and ^Di waves. 
Nevertheless, although there is a significant overall relative 
uncertainty between the tt" and tt^ production data, as- 
sociated with luminosity and other systematic effects, it is 
clear from the comparison of the seven free parameters in 
Eq. (5) with the ten observables in Table 1 that the system 
is over determined. As a consequence, if an acceptable solu- 
tion is achieved it would support the approximation made 
in our analysis, such as the neglect of higher partial waves, 
d-d interference and the effect of coupling between the 
and '^F2 partial waves 

The best fit to the combined pp — >■ {pp}s 7r° and 
pn — >■ {pp}s7r~ data sets is obtained with 

Mf= (55.3 ± 0.4) - (14.7 ± 0.1)i Vnb/sr, 

= -(26.6 ± 1.1) - (8.6 ± 0.4)z ^/nh/sI, 

Mf = 5.3 ±2.3 Vnb/sr, 

M^ = -(32.4 ± 2.2) + (17.3 ± 2.7)z \/nb/sr, 

M^ = -(109.6 ± 9.6) + (140.7 ± 4.0)i \/nb/sr. (6) 

Since this solution has x^/NDF = 89/82, it shows that our 
truncated expansions can give a very good description of 



the data. The contribution from the — >■ ^Sod transition 
is clearly very small and, if one eliminated M^ completely, 
it would give only a marginally poorer fit with x^/NDF = 
94/82. Note, whereas the phases of Mf , M[ and M[ were 
imposed (see above) , the phases of M^ and M^ were ex- 
tracted from the data. 

The quality of the fits can also be judged from the 
comparison of the curves in Figs. 3 and 4 with the data. 
The residual small discrepancies in the description of the 
analysing power might, of course, be due to the neglect of 
some of the smaller terms. However, it should also be borne 
in mind that the main systematic uncertainty, namely the 
relative normalisations between the pp — > {pp}s tt" and 
pn {pp}sTr^ data sets, has not been included in the de- 
termination of the parameters. On the other hand, adjust- 
ing this by a few percent would not lead to any qualitative 
changes in the solution of Eq. (6). 

Another way of judging the changes introduced by mak- 
ing a global fit to all the data simultaneously rather than fits 
to individual distributions is to look by how much the pa- 
rameters themselves were changed by this procedure. These 
are shown in Table 1. No changes at all are to be noticed 
for the pp — >■ {pp}sTT^ reaction and for the pn — {pp}s''^~ 
case it is only az{pn) where the difference is greater than 
the error bars. 

The conclusions that one can draw from Eq. (6) are first 
that, although d-wave pion production is significant, this 
is almost exclusively from the ^P2 state since the "^^2 — >■ 
^Sod transition is very weak. In the pn — )■ {pp}si^~ case the 
amplitudes are dominated by the ^Di — > ^Sop transition. 

A partial wave decomposition of the TRIUMF pn {pp}sT^ 
data was attempted [5] but at the time there were no 7r° 
production results available to constrain the fit and the 
authors had to rely on the application of the Watson the- 
orem even for the and ^Di waves, which are strongly 
coupled. Furthermore, their data did not extend over the 
whole angular domain and so it is not surprising that 
their partial wave results bear little relation to ours. They 
found negligible (i-wave production and did not identify 
the dominant role played by the ^Di ^S^p transition. 

In summary, we have measured the differential cross 
section and proton analysing power of 7t~ production in 
the quasi-free pn — >■ {pp}sTf~ reaction in the 353 MeV 
energy region over the whole angular domain. These com- 
plement the analogous data obtained on 7r° production in 
the pp — )• {pp}s TT*^ reaction. Through a careful use of the 
Watson theorem, an amplitude analysis of the combined 
data sets has been achieved and these have allowed the 
determination of the partial waves up to =2. 

Some of the biggest uncertainties in the global amplitude 
analysis arise from the normalisations in the pp — t- {pp}s 
and pn — >■ {pp}sT^~ data, which affect primarily the s- and 
p-wave production, respectively. The relative normalisation 
can be determined independently through a measurement 
of the transverse spin correlation parameter A^^x for the 
pn — >■ {pp}s7r~ reaction [21]. Such data have already been 
taken and, when the results are available in 2012, these will 



make the amplitude analysis more robust. 

The pion-production amplitudes extracted must now be 
analysed using chiral perturbation theory and this will be 
an important step to provide a deeper understanding of low 
energy pion dynamics. If this programme is successful, it 
will provide further strong evidence that chiral perturba- 
tion theory can indeed be applied to NN NNn provided 
that the large momentum transfers are taken into account 
properly. This is a precondition for the analysis of isospin- 
violating reactions in the same framework [22-24], which 
will shed light on the role of the quark mass term in low 
energy pion reactions. 
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